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We obtain a complete series solution of stationary Schrodinger's equation in the general quantum 
systems. It is exact in the sense that any approximation means is not used, or that the whole 
corrections or contributions from all order perturbations are involved if the perturbation concept is 
introduced. Furthermore, the useful forms of our exact solution and a new expression the complete 
Green operator are given out. As a universal and analytical solution, it is helpful for the theoretical 
^\ • derivations and practical calculations in quantum theory. 
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' Can we know an exact solution of stationary Schrodinger's equation (SSE) |l] if an arbitrary time-independent 
Hamiltonian in a general quantum system is given? The answer from the known theory of quantum mechanics is NO. 
But, our answer based on this article is YES. 

This problem is both basic in the theory of quantum mechanics and central in the applications of quantum mechanics. 
If the exact solution of SSE is known, then, the quantum state at any time is able to be exactly obtained for a given 
initial state, the different probabilities, the overlaps between the states and the expectation values of physical quantities 
. are able to be expressed exactly and calculated precisely, the non-perturbation frameworks of quantum mechanics 
i' is possible to be built, and some important conclusions and theorems should be able be derived out and proved. 
, However, the interesting on it may be getting a little weak because that this problem is too old, too difficult, and 
on the other hand, new developments and applications of quantum mechanics are very hot. At present, the exact 
solutions of SSE only exists in a few simpler quantum systems. For a general quantum system, this problem has 
not been solved since quantum mechanics emerged. In fact, this unsolved problem is an important origin of a lot of 
difficulties in theory developments and practical applications of quantum mechanics. It is unnecessary to say how 
important and how interesting of this problem because anyone working on quantum physics knows its answer. In this 
article, let us to study this problem in order to obtain an exact series solution of stationary Schrodinger's equation in 
the general quantum systems. 

Barring a few exceptions, the Schrodinger equation is always solved in a particular basis or representation. For a 
general quantum system, the form of its Hamiltonian H is not known. Thus, the chosen representation had better to 
■ be spanned by the set of orthogonal, normalization and complete eigenvectors of its solvable part, just as done in the 
perturbation theory and Green function method. That is, we split the Hamiltonian into two parts: 
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OV H = H Q + H 1 (1) 
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where the eigenvalues problem of Hq 

^ : #oi$ 7 ) = £ 7 |$ 7 ) (2) 
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is solvable. In other words, the eigenvectors |<I > 7 ) and the corresponding eigenvalues E 1 of Hq can be found or known, 
and the so-called Hq representation is constructed by the set {| ( i > 7 )}, where 7 takes over all possible indexes of energy 
levels. However, different from the usual perturbation theory, we do not introduce the perturbation concept. That is, it 
is unnecessary that Hi is small enough compared with Hq, or H\ is proportional to a very small perturbed parameter, 
because it will be seen that no perturbed parameter appears obviously in our solution, or all order corrections and 
contributions coming from H\ have been included for the case that H\ is indeed small enough compared with Hq. Of 
course, for this typical case, one also can take up some order approximation of perturbation in our following solution. 
However, the result is not trivially equivalent to the perturbative solution because that the reasonably dynamical 
rearrangement feature of our solution leads to the effects from the different cut approximation. 

It must be emphasized that our method to split the Hamiltonian of system into two parts has such a principle that 
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Hq matrix is diagonal and H\ matrix is off-diagonal in the representation spanned by {|<I> 7 }} , that is 

{%\H Q \^)=E^ rf , <«M#il*V> = {^ = (3) 

It is easy to prove that we always make such a division. If there are the non- vanishing diagonal elements of Hi matrix 
in the first division, we then redivide Hi matrix so that all of its non-vanishing diagonal elements are subtracted and 
then added to the Hq matrix in the first division, while all of its off-diagonal elements is taken as a new Hi matrix. 
It is obvious that the redivision result obeys our above principle. At this time, the eigenvalues of new Hq are added 
by the non- vanishing diagonal elements of the original Hi matrix, but the eigenvectors of new H are invariant. Note 
that we do not consider the case that all off-diagonal elements of Hi matrix are zero. It is trivial because that the 
eigenvalue problem of the Hamiltonian has been solved in this case. 

It should point out that an obvious reason that we split the total Hamiltonian matrix into diagonal part and off- 
diagonal part rather than the nonperturbed part and the perturbed part marked by a perturbed parameter in the 
normal perturbation theory is: we should, even have to, give up the perturbation concept in the beginning, otherwise 
we will stick in the difficulty same as the normal perturbation theory. Moreover, in our point of view, introducing 
the perturbed parameter and expanding according to its power is obviously a transcendent input, and this input covers 
the inherent feature that the perturbation series is multiple and excludes some important cases that the interaction 
strength is not very weak. In fact, giving up the perturbation concept is the most important reason why our exact 
series solution of SSE has new physics content. 

We first directly express our exact series solution of SSE just like an Ansatz so that it can be understood conveniently 
and easily. Then we prove it, and its derivation is arranged behind the verification. 

In a general quantum system, its stationary Schrodinger's equation reads 0] 

H\^ 1 )=E y \^ 1 ) (4) 
The exact series solution without any approximation of this equation is 

1*7) = l*7> + E Q^i') (5) 

where 

Q 7 , 7 = -J^ + E £ a .^-■-^ „ (6) 

E l - E f' 1=1 71 , 72 ,-, 7! ^7 i E 1 ~ E y)( E ~/ ~ E li)i E l ~ ^72) ' ' ' ( E l ~ E li) 

While the corresponding eigenvalue E~ as well as the quantity E 7 appearing in <3 7 < 7 definition is a real solution of 
following equation 

i? 7 (i? 7 — Ej) = Ej — E-y (7) 
Here, i? 7 (z) is called as a kernel function in our expression defined by 

R ( ) = V V g 77i g 7i72... g 7.7 

h 71 ,7 2 ^7 & - ^ - - E 12 - z) ■ ■ ■ (E, - Ell z) W 

Note that the convergence of our series solution is implied by the existence of real solution of Eq. . Otherwise we 
can not have above solution of SSE. 

Now let us directly prove that our solution is correct. Firstly, acting the total Hamiltonian on the right side of 
Eq.©, we get 

H (|* 7 )+ ]T Q 7 ' 7 |<IV} = £ 7 I*7>+E Q7'7^7'l^7')+ E5 7 ' 7 | $ V>+E E Q7'7.9 7 " 7 >7"> (9) 
\ 7'#7 / 7'#7 7'#7 7" 7'#7 

Secondly, we divide the fourth term into two parts with respect to 7" = 7 and 7" =/= 7 

E E ^'75 7 " 7 >7"> = Y,9 11 'Qi'i\^)+ E 5 7 " 7 'Q7'7l < V> (10) 

7" 7'#7 7'#7 7',7"#7 



Its first term can be rewritten as 
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By redenoting 7' by 71 in the first summation, and {7', 71, • • • ,7/} by {71,72, ■ • • ,7/+i} in the second summation, 
the above equation becomes 
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Then, we change the starting value of the summation in the second term of above equation via. the transformation 
I — > I + 1, we can arrive at 
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In its last step, Eq. ([7]) has been used. Similarly, using the trick of dummy index transformation, we have 
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= [E-y-E^jQ^-g-y 1 
Thus Eq. tfTO]) has the following form 



EE^ 7 " 7 >7"> = (X-£ 7 )i<i> 7 }+ E [( 



7" 7'#7 

Finally, substituting above relation into Eq.([9]), we obtain 
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It immediately goes to our conclusion. 

This proof is very easy even for a beginner in quantum mechanics if he/she is familiar to the method of dummy 
index transformation. Why such a complicated problem to solve the stationary Schrodinger's equation in a general 
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quantum system is changed to so simple? It is not because we can luckily guess this conclusion, but we obtain it 
based on our previous works [!, H[ • 

In the study on some applications of quantum mechanics, we often encounter even befall some difficult to understand 
and calculate problems, in particular, for open systems because the Hamiltonian is too complicated to solve its 
eigenvalue problem. Although the perturbative solution of SSE has been known, but its general term in c-number 
form has not been written obviously at our knowledge, only the Hamiltonian eigenvalues have a general term in the 
operator form. As to the Hamiltonian eigenvectors, their time-independent solution are expressed by Green function 
or operator in the potential scatter theory, as well as their dynamical form solution are written by the time evolution 
operator or the path integral. Consequently, we would like to seek a new way even formulism in order to solve the 
faced difficulties. 

As is understood by us, after giving up the perturbation concept, it is almost impossible to seek for the exact 
solution of the SSE directly. However, we know that the exact solution of dynamical Schrodinger equation can be 
obtained form the exact solution of SSE if the Hamiltonian is time-independent. Thus, vice versa, can we obtain the 
exact solution of SSE from the exact solution of dynamical Schrdinger equation? While the later in a general system 
has been known at least in form via the action of the time evolution operator. However, its obvious shortcoming is 
the expression of the time evolution operator is too formal to be usable directly in the theory and application. So, the 
key to seek for the exact solution of SSE is now changed to how to find out the general and explicit c-number form of 
time evolution operator. If this aim is arrived at, this new mentality of our research will have a successful possibility. 

As is well-known, to find the explicit c-number function form of the time evolution operator is a successful linchpin 
of Feynman path integral formulism and to expand the time evolution operator as a power series is a powerful 
headstream of Dyson series (in the interaction picture) [6]. Consequently, our physical idea [3] comes from the 
combination of Feynman path integral spirit and Dyson series kernel. We first derive out a general and explicit 
c-number expression of the time evolution operator that is different from them as the following: 

oo 

|*(t)> = £A™'(* y |*(0)>|* 7 >, A 11 ' = (<& 7 |e- ifft |<lv> = X>7 V (17) 

7,7' ;=o 

where the final state denoted by |^(t)), a given arbitrary initial state is denoted by | \T/ (0) ) and the coefficients AJ 1 
are defined by 

1+1 5 5 , ~-iE^t i 



^ =e ffi \', AT = £ E 7X 77 ' +1 lis™ 1 (18) 
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Note that it is still in a series form because that there is usually not the compact exact solution (or the exact solution 
in the closed form) of SSE for a general quantum system. It is clear that this expression can be thought of to be exact 
in the sense that any approximation means is not used, or that the whole contributions from all order perturbations 
are involved if the perturbation concept is introduced. 

It is very interesting to study the relation between our dynamical exact solution of the dynamical Schrodinger 
equation (JTTJ) and the results from the normal perturbation theory. By careful calculations and analyses, we found 
that our solution \YJ\ can naturally and reasonably involves partly corrections or contributions from the higher order 
perturbations [4j, but these corrections are cut out in the normal perturbation theory, for example, the transition 
probability in a very short time as well as the amplitude of nearer states in the perturbed wave function exist these 
problems [J,0|- Moreover, in the process of the calculations using the solution (TTT)) we observe two basic facts: (1) The 
total Hamiltonian should be split into diagonal and off-diagonal two parts rather than the perturbed and non-perturbed 
two parts; (2) Introducing perturbed parameter should be delayed as possible, had better do it in the final calculation. 
So we propose a Hamiltonian redivision skill and dynamical rearrangement technology in order to implement our above 
views. Furthermore, we build an improved scheme of perturbation theory based on our dynamical exact solution (JTTJ) 
so that the corrections or contributions from partly higher order perturbations are involved in the results, and so our 
method leads to that the related problems are more accurately and effectively calculated. Particularly, by analyzing 
our results, we propose a conjecture about the rearrangement series of Hamiltonian eigenvalues and realize a conclusion 
that the perturbed series of the Hamiltonian eigenvalues or Hamiltonian eigenvector is inherently multiple. In other 
words, the power expansion of perturbed parameter is not a best method. 

As an extension of our frameworks [4j, we continue to present the open system dynamics and the form of Green 
function and path integral J8|- In the study on them, we further understand why our results are needed. 

Obviously, there are many apparent singular points in the expression (|18|) of AJ 1 , but they are fake in fact. In our 
manuscript |4J this problem has been fixed by finding their limitations in terms of contraction and anti-contraction of 
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energy summation indexes. Recently, by theorizing this method and using partition function 
all the apparent singular points and arrive at 
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However, more importantly and interestingly, we obtain 
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It is proved by expending the partition function in Eq. (|19D into the time power series and verifying the coefficient 
power relation. In fact, the form of Eq. (|20[) has its physics origin rather than the mathematics arbitrariness, and it is 
valid in the general quantum systems independent of the form of Hamiltonian. Therefore, we think that the complete 
series expression of Hamiltonian eigenvalues in the general quantum systems is just: 
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Obviously, it is simply not a summation over the perturbed parameter, but a series of power of the kernel function 
R~/{z) as well as its derivative at z = 0. It is completely different from the normal perturbation theory in its 
thoughtway. In particular, when a cut is introduced in a practical calculation, a higher A 2 term than the last term 
of cut part is dropped, but not a higher A term than the last term of cut part is dropped in the normal perturbation 
theory. Moreover, by studying expression of Hamiltonian eigenvalues (HJ), we propose a calculating approach of 
eigenvalues of arbitrary Hamiltonian via solving an algebra equation (JT]) 3]. In addition, recalling the Richardson's 
exact solution in the pair model Q, we can see that there is a little similar in form between it and ours because the 
eigenvalues in the Richardson's exact solution also obeys an algebra equation, but the Rechardson's exact solution is 
only usable in a special system. 

Without loss of generality, we can expand the Hamiltonian eigenvector |\& 7 ) in the solvable representation spanned 



{|$ 7 }}, that is, 
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|$ 7 '). For simplicity, we set a 77 = 1 for a given 7 since that |^ 7 ) only can be 



determined to the difference from a normalization factor. Then, acting ($ 7 | on Eq.Q, we have E 7 = , g 11 a 77 '. 

It is easy to obtain its solution a 7 y = <5 77 < + Q 7 ' 7 (1 — ^77') if the total Hamiltonian eigenvalues is a solution of Eq.([7]). 
This derivation also verifies the result that the total Hamiltonian eigenvalues obey such an algebra equation ([7]) . 

Therefore, our solution of the stationary Schrodinger's equation in the general quantum systems is given originally 
based on our formal framework on quantum mechanics in the general quantum system, in particular, the conclusions in 
our paper Q , but it does not come from a guess although we can drop the details of our derivation. It can be thought 
of to be exact in the sense that any approximation means is not used, or that whole corrections and contributions from 
all order perturbations are involved if the perturbation concept is introduced. Breaking the accustomed mentality in 
the normal perturbation theory, finding the right division principle of Hamiltonian, using the more explicit c number 
expansion form of time evolution operator, dealing skillfully with and reasonably rearranging the perturbed series are 
advantages of our formulism f||. They bring us successfully to arrive at our purpose, and are also the reasons why 
our exact series solution of SSE can have more and new physics content. 

For convenience in the theoretical derivation and efficiency in the practical calculation, let us rewrite our solution 
in a more compact form. Firstly, introducing a project operator 



P 7 = /-|$ 7 )($ 7 | 

Here, / is an identity matrix and defining a limited operator of Hi 
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where we have used the facts that P 7 | ( i> 7 ) 



and ($ 7 |P 7 = 0. It implies that the elements of the 7-th row and 7-th 
column of i?i 7 matrix are zero, the other elements are same as the Hi matrix in the Hq representation. Then, we 
again define a revised Green operator 
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we arrive at 
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Note that Hi matrix is off-diagonal, that is, g 11 
identity is valid 



= (* 7 /|G 7 #i|$ 7 ) + ($ 7 '|G 7 #i]T (pyHi^ 1 G 1 H l \§ 1 ) (26) 
for any 7. If we focus on such quantum systems that the following 
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In fact, the same skill has been used extensively in the normal Green operator method, so it is still general enough. 
Therefore, the eigenvector of our exact series solution is rewritten as an operator form 



l*7> = 1*7) + E l^'XiM^ffi 

7' #7 



= \$> 1 )+P 1 G 1 H 1 
Comparing with the Green operator method 
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it immediately follows that the complete Green operator can be expressed by 
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When one tries to use it to the potential scatter problems, the boundary condition has to be considered. The related 
work is put in our another manuscript (in prepare). 

Similarly, the equation satisfied by the eigenvalues of our exact series solution can be rewritten as 



(^\Hi - G 7 gi|$ 7 ) = {^l^G = ~ #11*7) = E -r E i 



(32) 



This expression will largely advance the efficiency in the numerical calculations. Some examples studied by us are in 
progressing. 

It must be emphasized that our solution is a universal, systematical and programmable solution that is independent 
of the form of Hamiltonian. The preconditions of our solution is that the Hamiltonian contains a solvable part, which 
is the same as the perturbation theory and the Green function method. However, we need that the Hamiltonian is 
spilt into a diagonal Ho and a off-diagonal Hi, and it is always made. In fact, our solution can be thought of an 
Ansatz. Direct verification indicates its validity, and so its preconditions are looser. Of course, we have used the 
assumption of existence and uniqueness of solution of SSE for the general quantum systems. 

It is easy to verify that our results is consistent with ones in the normal perturbation theory if one expands our 
expressions of Hamiltonian eigenvectors and eigenvalues according to the order of perturbed parameter for the known 
lower order forms. It implies not only that our solution contains the perturbation solution, but also our solution 
have more physics content. As an exact series solution, our solution has a neat form of general term and involves the 
contributions from all order perturbation if the perturbation concept is needed. Actually, our solution is inherently 
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a non-perturbation solution. Comparing with the the Green function method, our solution does not need to solve 
the differential equation, and comparing with the directly diagonalized method, our expression is in a analytical 
form rather than a numerical form. Therefore, our solution is more suitable to the theoretical derivation and proof. 
Moreover, the convergence of our solution is considered in terms of the algebra equation satisfied by the Hamiltonian 
eigenvalues. In other words, the existence of solution of this equation implies that the series of our solution can 
converge to a finite number. 

A universal form has to pay the cost. A universal solution is unnecessary optimal for a special system. From the 
universality to speciality it still needs to study. In a concrete system, the form of our exact solution should be and can 
be simplified. In particular, we can not claim that our solution has the optimal efficiency in the numerical calculation 
for some special systems. This situation ought to be understandable, for example, the direct diagonalization method to 
find the eigenvectors and eigenvalues of Hamiltonian in some simple systems has the higher efficiency in the numerical 
calculations than the perturbed solution method to do this. However, ones cannot use this evidence to deny the 
importance of perturbative solution as well as perturbation theory. In fact, the direct diagonalization method is a 
good method limited within the numerical calculations. If there were not the perturbation theory and perturbative 
solution, will quantum mechanics have today achievement? The academic value and scientific significance of our exact 
series solution should be judged based on whether our solution contains more physics content than the perturbative 
solution. 

Apparently, our solution seems a rearranging summation when the perturbation concept is introduced, but essen- 
tially it reveals the physics nature in a quantum system. This is an important reason why a complicated problem 
to solve the stationary Schrodinger equation in a general quantum system is changed to so simple. In fact, our 
solution directly comes from quantum dynamics and it has the completeness, orderliness and clearness. In practical 
calculations, they also provide a physics method how to choose part corrections or contributions from higher order 
perturbations, which is able to simplify the calculation and lead the result more precise. Therefore, our solution is not 
trivially equivalent to the normal perturbative solution when the perturbation concept is introduced because that the 
reasonably dynamical rearrangement feature of our solution leads to the effects from the different cut approximation. 
Two examples have been presented in Refs.[H, 0]- In special, using the expressions of our exact solution of SSE 
and (|32[) for many interesting systems, which satisfies the condition (|27p . we can obtain some encouraging results. 

By comparison and analyses, our results are helpful for theoretical derivations and practical calculations since its 
universal, analytical and exact features. However, it must be pointed out that we expect that our exact series solution 
of SSE has more important and interesting conclusions and applications. It will be more clearly seen the academic 
value and scientific significance of our exact series solution of SSE. 

For simplicity, we only consider the non-degeneracy and discrete case here, but our derivation can be extended to 
the degeneracy and/or continuous case. At present, we are studying some concrete examples in order to further verify 
the validity of our results and reveal the ability of our approach. We are sure that it will bring many interesting 
conclusions when our frameworks [H, 3, 0> @] & r e extended to the time-dependent cases and quantum field theory. 

I am very grateful for Zhou Li for her irreplaceable contributions from the expression and calculation of Hamiltonian 
eigenvalues in the general quantum systems [|[ . This work was supported by the National Natural Science Foundation 
of China under Grant No. 10975125. 
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